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ABSTRACT 

We present the first extensive photometric results of CL Aur from our BVRI CCD 
photometry made on 22 nights from 2003 November through 2005 February. Fifteen 
new timings of minimum hght were obtained. During the past 104 years, the orbital pe- 
riod has varied due to a periodic oscillation superposed on a continuous period increase. 
The period and semi-amplitude of the oscillation are about 21.6 yrs and 0.0133 d, re- 
spectively. This detail is interpreted as a light-travel-time effect due to a low-luminosity 
K-type star gravitationally bound to the CL Aur close system. Our photometric study 
indicates that CL Aur is a relatively short-period Algol-type binary with values of 
g=0.602 and i=88°.2. Mass transfer from the secondary to the primary eclipsing com- 
ponent is at least partly responsible for the observed secular period change with a rate 
of dP/dt = +1.4xl0~^ d yr~^. A cool spot model has been calculated but we think 
that an alternative hot-spot model resulting from a gas stream impact on the hot star 
is more reasonable despite two difficulties with the explanation. Absolute dimensions 
of the eclipsing system are deduced and its present state is compared with tracks for 
single star and conservative close binary evolution. Finally, we examine the possible 
reconciliation of two different calculations of the luminosity of the hot spot and a re- 
interpretation of the secular term of the period variability. 



Subject headings: binaries: close — binaries: eclipsing — stars: individual (CL Aurigae) 
— stars: spots 
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1. INTRODUCTION 

Algol-type close binaries are semi-detached interacting systems in which one type of interaction 
is mass transfer between the component stars by means of a gas stream. They have been known 
as good astrophysical laboratories for studying accretion processes because a number of them are 
bright. They are in the slow phase of mass transfer with dM/dt ~ 10^^"*^ — 10^'^ Mq yr~^ and 
do not undergo violent eruptions that interfere with the accretion process. The circumstellar 
structures produced by the mass-transfer process in these systems have been sorted according to 
orbital period by Richards &: Albright (1999) but do not depend upon it significantly. Rather, their 
natures can be easily understood from the position of the mass-gaining component in the so-called 
r-q diagram in which the fractional radius r = {R/a) of a gainer is plotted versus the mass ratio 
q and compared with the semianalytical computations of the gas stream hydrodynamics of Lubow 
& Shu (1975). In the short-period Algols located above the Ud curve of the diagram (cf. Figure 2 
of Richards & Albright), the hot, detached primary star is large relative to the orbital radius and 
the two components are too close to each other to form an accretion disk or even a stable accretion 
annulus. Instead, it is possible that an impact region, and hence a hot spot, can be formed on 
the surface of the primary star somewhat displaced from the line of centers due to the Coriolis 
acceleration imposed on the flowing gas. If the secondary stars are sufficiently cool, they likely 
display enhanced magnetic activity due to deep outer convective layers and rapid rotation. This 
magnetic mechanism may contribute to the period and light variations for systems with spectra 
later than F-type (Hah 1989). 

CL Aur (GSC 2393-1455, HV 6886, TYC 2393-1455-1) was discovered to be a variable star by 
Hoffleit (1935) based on photographic plate estimates. Kurochkin (1951) presented the first (partial) 
photographic light curve of the star and the original light elements, Min. I = HJD 2,432,967.262 -|- 
1.2443666E'. The value of the period positions this object toward the short-period limit for Algols. 
The spectral type of the primary star was classified to be AO by Gotz &; Wenzel (1968). Since 
then, times of minimum light have been published assiduously by numerous workers but, to our 
knowledge, a complete light curve and the fundamental parameters for the binary system have not 
been made so far. Changes of the orbital period have been considered by Hegediis (1988) and Wolf 
et al. (1999). Hegediis selected this system as a possible candidate for the study of apsidal motion. 
However, the later authors ruled out this possibility from CCD timings for primary and secondary 
eclipses. They suggested the cause of period variation to be a light-travel-time (LTT) effect due 
to the presence of a third body in the binary system. Most recently, Wolf et al. (2007, hereafter 
W07) reported that a long-term period increase is superimposed on an LTT orbit with a period of 
P3=21.7 yrs, a semi-amplitude of i^=0.014 d, and an eccentricity of e=0.32. 

In the Simbad data bas^, the system is described as an eclipsing binary of (3 Lyr type. 
BVJHK magnitues are listed for the star but these are from heterogeneous sources and are not mu- 



^http://simbad. u-strasbg.fr/simbad/ 



-3- 



tually consistent. Part of this inconsistency arises because the magnitudes refer to different phases 
in the Keplerian cycle. With the weh-known transformations (ESA 1997), standard photometric 
vahics for CL Aur in the Johnson system were calculated to be y = +11.62 and (B — V) = +0.33 
from the Tycho results. These refer to some unknown Keplerian phase. These are not consistent 
with those in Simbad presumably because of the large phase-locked variations in magnitude and 
color index of the binary. 

At present, CL Aur is known only as a neglected eclipsing system composed of an A- type 
primary and a cooler companion. In order to derive photometric solutions and to examine whether 
the WOT suggestion is appropriate for the orbital period change, we decided to obtain light curves 
with multiband photometry. In this paper, we present the first mutual analyses of the O-C diagram 
and the light curves. 

2. CCD PHOTOMETRIC OBSERVATIONS 

New photometric observations of CL Aur were obtained using a SITc 2K CCD camera and 
a BVRI filter set attached to the 61-cm reflector at Sobaeksan Optical Astronomy Observatory 
(SOAO) in Korea. The observations of the first season were made on 14 nights from 2003 November 
to 2004 March and those of the second season on 8 nights from 2004 December through 2005 
February. The exposure times were 75—140 s for B, 45—85 s for F, 33—65 s for R, and 30—60 s 
for I, respectively, depending on weather conditions. The instrument and reduction method have 
been described by Lcc ct al. (2007) and a 2x2 binning mode was selected. The nearby stars GSC 
2393-1424 and GSC 2393-1418, imaged on the chip at the same time as the variable, were selected 
as comparison and check stars, respectively. Coordinates and Tycho magnitudes for the three stars 
of interest are given in Table 1. The differential atmospheric extinction among the three stars is 
negligible within observational error. Measurements of the check and comparison stars indicate that 
the latter remained constant throughout the observing interval. The Icr-values of the dispersions 
of the magnitude differences between them are about ±0.01 mag for all bandpasses. 

A total of 2747 individual observations was obtained among the four bandpasses (711 in B, 693 
in y, 685 in R, and 658 in /) and a sample of them is listed in Table 2. The light curves of CL Aur 
defined by our CCD photometry are plotted in Figure 1 as differential magnitudes versus orbital 
phase, which was computed according to the ephemeris for our hot-spot model determined later in 
this article with the Wilson-Devinney synthesis code (Wilson k, Devinney 1971, hereafter W-D). 
The filled and open circles are the individual measures of the first and second observing seasons, 
respectively. Our light curves from the first observing season had not defined secondary eclipse 
adequately. Except formally for the B bandpass, mean brightness differences between the two 
seasons (in the sense of Season 1 minus Season 2) are constant and smaller than the observational 
error of ±0.01 mag: +0.013±0.014 mag for B, +0.007±0.016 mag for V, +0.001±0.017 mag for R, 
and — 0.009±0.015 mag for /, respectively. Although we show no figure illustrating variability of the 
color indices, their phase-locked variations are large and in the expected senses. These variations 
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convey the idea of eclipses that are complete or very nearly so. 

In addition, we observed several eclipse timings using a SBIG ST-8 CCD camera attached to the 
35-cm reflector at the campus station of the Chungbuk National University Observatory (CbNUO) 
in Korea. The observations were made without a filter and reduced with the conventional IRAF 
package. The details of the CbNUO observations have been given by Kim et al. (2006). 

3. ORBITAL PERIOD STUDY 

We determined fifteen times of minimum light from all our CCD observations using the method 
of Kwee &: van Woerden (1956, hereafter KvW). These timings are listed in Table 3 together with 
all other CCD timings; the SOAO data are weighted means from the observations in the BVRI 
bandpasses. The second column gives the standard deviation of each timing. For the period study 
of CL Aur, 198 archival timings (16 plate, 101 visual, 8 photographic and 73 CCD) have been 
collected from the literature and from our measurements. Most of the earlier timings were actually 
extracted from the data base published by Kreiner et al. (2001). Because almost all but the CCD 
timings were published without error information, the following standard deviations were assigned 
to timing residuals based on observational method: ±0.022 d for sky-patrol plate or photographic 
minima, and ±0.007 d for visual minima. Relative weights were then calculated as the inverse 
squares of these values consistent with the errors and weights for the CCD timings. 

In order to examine whether the period change of CL Aur can be produced by a quadratic 
plus LTT ephemeris as suggested by W07, we fitted all times of minimum light to that ephemeris 
form: 

C = To + PE + AE^ + T3, (1) 

where r3 symbolizes the LTT effect due to a third body (Irwin 1952, 1959) and includes five 
parameters (ai2sini3, e, u, n and T). Here, n and T denote the Keplerian mean motion of the 
mass center of the eclipsing pair and the epoch of its periastron passage, respectively. In this 
analysis, the Levenberg-Marquart technique (Press et al. 1992) was applied to solve for the eight 
unknown parameters of the ephemeris. By using the third-body parameters of W07 as initial values, 
we obtained improved (in the sense of errors smaller than formerly) parameters for them and the 
results are listed in Table 4 together with other related quantities. In this table, P3 and K indicate 
the cycle length and semi-amplitude of the LTT orbit, respectively, and /(M3) is the mass function 
of the system. Within errors, our LTT parameters are not significantly different from those of W07. 
The sample masses (M3 sinis) of the assumed third body in the table were calculated by using the 
absolute dimensions of the eclipsing pair presented in a later section. 

The O-C diagram of CL Aur constructed with the linear terms in Table 4 is drawn in the 
top panel of Figure 2. The timings are marked by different symbols according to observational 
method and type of eclipse. The continuous curve and the dashed, parabolic one represent the 
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full contribution and the quadratic term of the equation, respectively. The middle panel displays 
the residuals from the linear and quadratic terms, and the bottom panel the CCD residuals from 
the full ephemeris. These appear as O-Cfuu in the fourth column of Table 3. As indicated by the 
figure, the quadratic plus LTT ephemeris gives a satisfactory representation of the ensemble of the 
residuals. In addition, another modulaton with a period of about five yrs and a semi-amplitude 
of about 0.0004 d may exist in the O-Cfuii residuals. A large number of future accurate timings 
is required before this can be tested at an acceptable level. If the orbit of the outer component 
is coplanar with that of the close eclipsing pair (13=88°. 2), its mass is about M3=0.83 Mq and 
corresponds to a spectral type of K1-K2 for a normal main-sequence star. This would contribute 
about 0.9% to the total light of the triple system, so it will be difficult to detect such a third light 
source from light-curve analysis. 

The quadratic term (A) of the ephemeris indicates a continuous period increase at a rate of 
-|-(1.4±0.2) X 10~^ d yr"^. Because CL Aur is a semi-detached system with the less massive and cool 
secondary component filling its inner Roche lobe (cf. Section 4), its Roche-geometry configuration 
permits mass transfer from the secondary to the more massive primary star, and it is conventional 

to ascribe such an increase to conservative mass transfer between the stars in the system. The 
calculated mass transfer rate is of the order of 1.3x10"^ Mq yr~^, among the largest rates for 
Algol-type systems. 

Alternatively, the 21.6 yr oscillation in the 0-C residuals could be caused by magnetic mod- 
ulation due to an activity cycle in the convective envelope of the late-type star (Applegate 1992, 
Lanza et al. 1998). The hot primary component of CL Aur likely has a radiative envelope as 
surmised from its spectral type while the less massive and cool secondary should have a shallow 
convective shell and at most weak magnetic activity. We applied the period (P3) and amplitude 
(K) to Applegate's formulae and obtained model parameters for possible magnetic activity. The 
results are listed in Table 5, where the bolometric magnitude difference (Amj-ms) relative to the 
mean light level was obtained with equation (4) in the paper of Kim et al. (1997). Most of the tabu- 
lated values are close to those derived for several other binaries that appear to support Applegate's 
theory, but the light variation predicted from an active CL Aur secondary is at the upper limit of 
the theoretical value (AL/L2 ~ 0.1) proposed by him. Since our observations represent the first 
complete light curve, we cannot test whether such a variation has occurred in the past. However, 
because the secondary star is expected not to be strongly magnetically active, we think the most 
probable explanation of the periodic oscillation to be the LTT effect due to a low-luminosity K-type 
tertiary companion. 



4. LIGHT-CURVE SYNTHESIS 

As shown in Figure 1, our observations clearly indicate that the light curve morphology of CL 
Aur is not jS Lyr type but rather very similar to that of Algol, its class prototype. To understand 
the geometrical structure and the physical parameters of the system, our BVRI light curves were 
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solved simultaneously in a manner similar to those for XX Cep (Lee et al. 2007) and GW Gem (Lee 
et al. 2009a) by using the latest versior^l of the W-D code and an extensive g-search procedure. 
The surface temperature of the primary star was fixed at Ti =9,420 K, according to its spectral 
type AO and Harmanec's (1988) table. We had attempted to improve the spectral classification by 
obtaining high-resolution spectra with the echelle spectrograph attached to the 1.8-m telescope at 
Bohyunsan Optical Astronomical Observatory (BOAO) in Korea (Kim et al. 2002). These images 
were not well enough exposed to improve upon the literature value but indicated a similar result. 
Initial bolometric and monochromatic limb-darkening coefficients were taken from the tables of 
van Hamme (1993) and were used together with the model atmosphere option. The (^-search for 
modes 2, 3, 4 and 5 of the synthesizing code (Wilson & Biermann 1976) converged and showed 
acceptable photometric solutions only for mode 5 (semi-detached systems for which the secondary 
stars fill the inner Roche lobes). As displayed in Figure 3, the optimal solution is around g=0.60. 
This undertanding conforms to the sense of the mass transfer from the secondary component to 
the hotter, more massive primary star suggested by the period study. 

The q value was treated as an adjustable parameter in all subsequent syntheses deriving binary 
parameters. The best result for unspotted photospheres is listed in columns (2) and (3) of Table 6 
and the residuals from the analysis are plotted in the left panels of Figure 4, where phase-locked, 
unmodelled light variations are indicated. Such features can be reasonably attributed to a hot 
spot on the surface of the primary as a result of the impact of the gas stream from the cooler, less 
massive secondary star. This situation is known for the Algol-type, semi-detached binary RZ Cas 
(Rodriguez et al. 2004). 

This interpretation does not exclude the possible existence of magnetic cool spots located on 
the surface of the late-type star. We therefore tested two spot models: a hot spot on the more 
massive primary star due to mass transfer and a cool spot on the secondary star caused by magnetic 
activity. In a formal sense, as shown by the entries on the last line of Table 6, the hot spot model 
does improve the light-curve fit. Separate trials for a cool spot on the secondary star were not so 
successful as for the hot-spot model. Final results for all the light curves are given in columns (4) 
and (5) of Table 6 and the residuals from our hot-spot model are shown in the right panels of Figure 
4. As seen in the figure, there is a slight improvement in the residuals for the spot model at about 
phase 0.87 in all panels compared to the unspotted model. This orbital phase is almost exactly 
where we would expect the effects of a gas stream to be evident. The hot spot agrees well with the 
concept of mass transfer from the secondary to the primary component inferred from the period 
analysis and from the relatively large size of the hot star. In addition, there are small systematic 
differences between the two seasons for the light residuals, which means that mass transfer activity 
may be variable, as is sometimes reported for light curves of Algol binaries. In all these trials, a 
possible third light source (£3) was considered but the results remained indistinguishable from 0.00 
within their errors. We fixed £3 to be 0.0 during the final light-curve analysis. 



^ftp: //ftp. astro. ufi.edu/pub/wilson/ 
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Because of the almost complete eclipses, the light curve determinacy for CL Aur is quite high. 

In idealized Roche geometry, the impact of the gas stream should be on the primary star's equator 
but this is not the case for CL Aur. We return to this detail later. It could also be true that both 
hot and cool spots exist and that our data arc sufficient to isolate only the dominant one. 

A recent study by Lcc ct al. (2009b) indicates that the minimum epochs of the cool contact 
binary AR Boo have been systematically shifted by light asymmetries due to spot activity, as 
predicted by Maceroni k, van't Veer (1994). To check this possibility for this Algol-type system, we 
calculated the timings for each of our CL Aur eclipses with the W-D code. The results are listed in 
the second column of Table 7, together with the minimum times obtained by the KvW method for 
comparison. We can see that the differences among them in column (4) of Table 7 are within the 
precision of each KvW minimum and the hot spot does not inflect those timings. This agreement is 
doubtless due to the small temperature contrast of the spot against the photosphere and the small 
angular deviation of the spot from the systemic line of centers. 

A new class of 6 Set stars designated as oEA (oscillating EA) objects (Mkrtichian et al. 
2004) has been identified as the (B)A-F spectral type mass-gaining components of Algol-type semi- 
detached systems. About half of them have been discovered through a photometric survey project 
run by a Korean group (cf. Kim et al. 2003) in order to search for A-type pulsating components 
in classical Algols (Soydugan et al. 2006). The oEA stars have pulsational amplitudes and periods 
similar to classical 5 Set stars but a different evolutionary scenario due to mass accretion. From 
its spectral type, the primary star of CL Aur would be a candidate for such pulsations. Therefore, 
we applied the discrete fourier transform program PERIOD04 (Lenz &; Breger 2005) to the light 
curve residuals from our binary models, but there was found no periodicity with a semi-amplitude 
larger than 3 mmag. 

5. ESTIMATED ABSOLUTE DIMENSIONS 

The foregoing presents a consistent picture of CL Aur in the sense that a determinate repre- 
sentation of the light curve has been achieved and that the secular period variability has a highly 
probable cause consistent with the light curve interpretation. Now it is necessary to look at further 
consequences of this description. 

Absolute dimensions of CL Aur can be estimated from our photometric solutions with the 
hot-spot model in Table 6 and Harmanec's relation between spectral type and mass. By assuming 
the primary star to be a normal main-sequence one with a spectral type of AO V, the astrophysical 
parameters for the components were obtained to be those listed in Table 8, where the radii are 
the mean volume radii calculated from the tables of Mochnacki (1984). The luminosities and 
the bolometric magnitudes were computed by adopting reff0=5,78O K and Mboi0=+4.69. The 
bolometric corrections were obtained from the relation between log T and BC given by Kang et al. 
(2007). The intrinsic color of {B — V)q = +0.06 for the binary system was estimated from their 
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calibration {B — V)i versus Ti of the primary star and from our computed light ratio; it leads to 
E{B — V)=+0.28. For lack of better evidence, it is prudent to assume that the mean system color 
index and interstellar reddening can be in error by as much as ±0.05. With this qualification and 
with the values of V and Mv,totai and the interstellar extinction of Ay=3.1E{B — V), we have 
calculated a nominal distance to the system of about 1,150 pc and an accidental error associated 
this distance of about 5 %. This result is also misleading for another reason: one does not know the 
phase at which the Tycho measures were made so there is no assurance that they refer to maximum 
light. The total error on the distance determination may well be 20%. 

From the estimated parameters of CL Aur, it is possible to consider the evolutionary state of 
the eclipsing system in mass-radius and mass-luminosity diagrams. The primary star lies in the 
main-sequence band between the ZAMS and the TAMS, while the secondary is larger and brighter 
than expected for its mass. The locations of the components in the Hertzsprung-Russell (HR) 
diagram are shown in Figure 5, together with single-star evolutionary tracks having masses of 2.20 
Mq and 1.35 Mq and solar metallicity (Girardi et al. 2000). There also appear in the figure 
theoretical mass-conserving evolutionary tracks for model donor (3.0 Mq, the present secondary) 
and gainer (1.2 Mq) stars belonging to a system with a total mass of 4.2 Mq and an initial orbital 
period of 1.0 d calculated by the Brussels groupH (De Loore & van Rensbergen, 2005). Although 
these tracks are not for a system with the same mass ratio and total mass as CL Aur (3.6 Mq), 
the two components can be scaled to fit the concept of a similar conservative model quite well, 
indicating that the binary system may still be undergoing mass transfer as a result of case A 
evolution. The present age of the system is estimated to be about 0.24 Gyr. 



6. DISCUSSION AND CONCLUSIONS 

We next probed the credibility of the hot (inferentially an impact) spot model motivated largely 
by two recognitions: (a) the two arrays of light curve residuals in Figure 4 do not appear very 
different one from another although formally the light curve fitting is improved by the assumption 
of a spot and (b) the center of the spot is sensibly distant from the stellar equator and orbital 
plane. 

The first matter is the consequence of the envelope of the residuals responding weakly to the 
spot modeling although the residuals of small absolute value did respond well with considerably 
reduced values on average. This can only be construed as a situation in which there exist residuals 
of noise much greater than most of the rest of them. In principle, these residuals could be produced 
by a combination of physical causes such as magnetic activity from the cool secondary and a change 
in the mass transfer rate. However, the actual noise level is not significantly larger than that (±0.01 
mag) of SOAO data made during the last few observing seasons and probably can be traced to 
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marginal observing conditions on some nights. 

The location of the spot has more physical interest. First, an extensive set of restricted 3- 
body calculations examined the longitude coordinate of the spot on the hot primary for a range 
(0.0000003 km - 450 km s~^) of initial velocities streaming from the Li point as shown in the 
first column of Table 9. This range obviously includes the thermal velocity in the cool star envelope 
with the lower limit of the range essentially that of free-fall. For the greatest extent of this velocity 
range the impact spot is close (within 21°) to the line of centers. In the hydrodynamic computations 
of Lubow k, Shu (1975), just such an effect appears with the streaming gas deflected about 20° 
from the line of centers, not nearly enough to avoid contact with the large, detached primary. 
From this point of view, the credibility of the concept of impact is therefore high, and conservative 
mass transfer of the streaming material must be considered likely. These calculations also make it 
difficult to imagine any stable accretion structures around the hot star but they invariably led to 
an orbital-plane spot rather than the modeled one of co-latitude very different from 90°. Another 
suite of calculations, also given in Table 9, showed the expected result: inital non-zero-velocity 
z-componcnts in the streaming gas led to impacts away from the hot stellar equator and these did 
not have to be large in order to fall 15-20° away from the equator. Is this result to be taken as 
evidence that the spot really exists at its modeled location or is it just a dynamical truism and 
doesn't verify the spot existence at all? The resolution of this quandary requires a mechanism 
to move the gas appropriately and three possible ones come to mind, (a) Turbulence in the gas 
moving from the Li-point caused a concentration of it to fall at the modeled position during the two 
observing seasons, (b) A weak magnetic field seated in the cool star had at least one component 
channeling the ionized fraction of the streaming gas to the modeled impact spot, (c) The solution 
for the spot is possibly not so unique as W-D indicates. For want of independent evidence, each of 
these is an unprovable hypothesis. 

More information does exist, however. We calculated the impact luminosity from the stream 
at the rate of mass transfer given by the interpretation of the 0—C diagram. For the same range 
(and in the same sense as above) of stream velocity, the impact luminosity varied between 3.2 Lq 
and 0.7 Lq. If the impact energy is partitioned between virialization into the gainer star and spot 
luminosity as Hilditch (1989) proposes, the least luminous output from a spot should therefore be 
about 0.35 Lq. This is to be compared to the spot luminosity from its black body W-D model 
of about 4xl0~^ Lq. The magnitude of the discrepancy means cither (a) that conversion from 
kinetic to luminous energy is very inefficient or (b) that the mass-motions rate indicated by the 
O-C diagram is not restricted to material leaving the Li point or (c) some combination of these 
two ideas or (d) that the concept of a spot is itself incorrect. For (a) to be a realistic interpretation, 
deep penetration of the impacting gas into the hot star is required so that most of the impacting gas 
becomes thermalized in the hot star. It is also necessary to postulate an envelope circulation pattern 
that re-surfaces at the modeled spot position for the residual gas that has not been thermalized. 
If (b) is to be entertained seriously, there must be substantial mass lost by evaporation from the 
entire Roche lobe that is the photosphere of the cool star. This demands only a small (~ 10~^) 
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contribution to the O-C diagram from the gas leaving Li with the majority of the mass and angular 
momentum loss contributed by that departing the rest of the cool star's surface. On a small scale 
the transfer would be conservative as usual but large scale mass loss would be the dominant cause 
of the ^-term in the ephemeris. Hypothesis (d) would be the most economical interpretation but 
conflicts with the impersonal syntheses of the light curves. At this time, we favor a combination of 
interpretations (a) and (b). 

In summary, our study of the orbital period and the light curves reveals that CL Aur is a 
classical Algol-type interacting system with the less massive and cool secondary star filling its 
inner Roche lobe. The possibility of a hot-spot model due to impact of streaming gas onto the 
hot star has led to confusing difficulties that we have resolved only tentatively. High-resolution, 
near-IR spectroscopy should reveal the lines of the secondary star and lead to accurate absolute 
parameters to replace our estimates. Moreover, there is also the possibility of obtaining direct 
evidence of mass-transfer activity, such as complex and variable line profiles of Hcc or variations in 
the strength of the O I 7774 absorption line. The evolutionary status of the system will then be 
more convincingly in hand. 

The more general reality of things is that there are many short-period Algols which have not 
been studied so deeply as this present work has examined CL Aur. It would make a significant 
advance if more of these - having different algebraic signs for the secular term of the period vari- 
ability - could be brought to the same level of knowledge as the present binary. For instance, should 
all such binaries with ^ > require a spot on the hot star near the line of centers, there would 
be major support for the reasoning concerning mass movements. Should it also happen that the 
spot luminosity was consistently found to be lower than the kinetic impact required for 50% energy 
conversion, there would be good reason to believe that the majority of the mass lost from the cool 
secondaries is, in fact, lost to the systems and not transferred conservatively. The short-period 
binaries found to require mode 4 (semi-detached systems with the primary stars filling its inner 
Roche lobes) for their representations would be expected, then, to show small-scale mass transfer 
to the secondaries and systemic mass loss from the primaries. Much valuable observational and 
modeling work remains. 
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Fig. 1. — SOAO light curves of CL Aur in the B, V, R, and I bandpasses. The filled and open 
circles are the individual measures of the first (2003 Nov— 2004 Mar) and second (2004 Dec— 2005 
Feb) seasons, respectively. Because of the high density of the points, many of the 2004 measures 
cannot be seen individually. 
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Fig. 2. — The 0-C diagram of CL Aur constructed with the hnear terms in Table 4. In the top 
panel, the continuous curve and the dashed, parabolic one represent the full contribution and the 
quadratic term of the equation, respectively. The middle panel displays the residuals from the linear 
and quadratic terms, and the bottom panel the CCD residuals from the complete ephemeris. CC, 
VI, PG, and P denote CCD, visual, photographic, and photographic plate minima, respectively. 




Fig. 3. — The behavior of S (the sum of the residuals squared) of CL Aur as a function of mass 
ratio showing a minimum value near g=0.60 for Mode 5 of the W-D code. 
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Fig. 4. — The magnitude residuals in B, V, R and / corresponding to two binary models in 
columns (2)-(5) of Table 6: without (left panels) and with (right panels) a hot spot on the primary 
component. The symbols are plotted in the same sense as Figure 1. 
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Fig. 5. — The HR diagram of the primary (P) and secondary (S) stars of CL Aur. The dashed 
hnes correspond to single star evolutionary tracks for the present masses of both components. The 
green- and pink-colored lines represent theoretical tracks for closely similar donor and gainer stars, 
respectively. See the text for details. 
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Table 1. Coordinates and photometric data for the program stars. 



Star 


GSC 


RA (J2000) 


DEC (J2000) 


V* 


{B - F)^ 


CL Aurigae 


2393-1455 


05'^12'^54!19 


-F33°30'28'.'4 


-hll.65 


+0.39 


Comparison 


2393-1424 


05'^13'^27!48 


-F33°26'46'.'3 






Check 


2393-1418 


05'^12"^19?08 


-h33°26'31'.'5 


-hl2.30 


+0.14 



tProm the Tycho-2 Catalogue (H0g et al. 2000). 



Table 2. SOAO CCD photometric observations of CL Aur. 



HJD 


AB 


HJD 


AV 


HJD 


AR 


HJD 


AI 


2,452,974.23806 


0.707 


2,452,974.23966 


0.525 


2,452,974.24094 


0.348 


2,452,974.23583 


0.104 


2,452,974.24426 


0.872 


2,452,974.24587 


0.673 


2,452,974.24716 


0.452 


2,452,974.24215 


0.179 


2,452,974.24997 


1.020 


2,452,974.25159 


0.775 


2,452,974.25288 


0.532 


2,452,974.24832 


0.280 


2,452,974.25568 


1.042 


2,452,974.25732 


0.760 


2,452,974.25861 


0.496 


2,452,974.25404 


0.312 


2,452,974.26142 


0.949 


2,452,974.26304 


0.682 


2,452,974.26433 


0.412 


2,452,974.25977 


0.281 


2,452,974.26721 


0.800 


2,452,974.26882 


0.540 


2,452,974.27010 


0.305 


2,452,974.26550 


0.196 


2,452,974.27288 


0.646 


2,452,974.27450 


0.420 


2,452,974.27579 


0.193 


2,452,974.27127 


0.106 


2,452,974.27860 


0.507 


2,452,974.28020 


0.299 


2,452,974.28149 


0.098 


2,452,974.27697 


0.030 


2,452,974.28423 


0.377 


2,452,974.28585 


0.198 


2,452,974.28714 


0.014 


2,452,974.28266 


-0.049 


2,452,974.29000 


0.255 


2,452,974.29161 


0.088 


2,452,974.29289 


-0.082 


2,452,974.28831 


-0.123 



Note. — This table is available in its entirety in machine-readable and Virtual Observatory (VO) forms in 
the online journal and also at the Web page ( http: / / binary. cbnu.ac.kr /bbs/zboard.php?id=lab_photometry| . 
A portion is shown here for guidance regarding its form and content. 
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Table 3. Observed CCD times of minimum light for CL Aur. 



HJD 




Error 


Epoch 


O-Cfuu Min 


References 


2,450,044 


3958 


±0 


0007 


-42 


5 


-0 


00122 II 


Wolf et al. (1999) 


2,450,097 


2826 


±0 


0002 








-0 


00014 I 


Wolf et al. (1999) 


2,450,714 


4888 


±0 


0006 


496 








00020 I 


Wolf et al. (1999) 


2,450,831 


4589 


±0 


0001 


590 








00019 I 


Wolf et al. (1999) 


2,450,884 


3431 


±0 


0003 


632 


5 


-0 


00100 II 


Wolf et al. (1999) 


2,451,157 


4815 


±0 


0006 


852 








00027 I 


Wolf et al. (1999) 


2,451,177 


3914 


±0 


0001 


868 








00043 I 


Wolf et al. (1999) 


2,451,508 


391 


±0 


0030 


1134 








00080 I 


Blattler (2000) 


2,451,570 


6077 


±0 


0003 


1184 





-0 


00037 I 


Baldwin & Samolyk(2004) 


2,451,880 


4526 


±0 


0002 


1433 





-0 


00061 I 


Brat et al. (2009) 


2,451,901 


6067 


±0 


0002 


1450 





-0 


00062 I 


Wolf et al. (2007) 


2,451,901 


6075 


±0 


0003 


1450 








00018 I 


Brat et al. (2009) 


2,451,921 


5164 


±0 


0002 


1466 





-0 


00068 I 


Agerer & Hiibscher (2002) 


2,451,925 


2501 


±0 


0021 


1469 





-0 


00006 I 


Brat et al. (2007) 


2,451,956 


3593 


±0 


0015 


1494 








00014 I 


Brat et al. (2007) 


2,452,013 


5998 


±0 


0002 


1540 








00004 I 


Baldwin & Samolyk(2004) 


2,452,017 


3345 


±0 


0003 


1543 








00165 I 


Wolf et al. (2007) 


2,452,196 


5200 


±0 


0003 


1687 





-0 


00113 I 


Diethelm (2001) 


2,452,252 


5171 


±0 


0001 


1732 





-0 


00047 I 


Wolf et al. (2007) 


2,452,333 


4014 


±0 


0001 


1797 





-0 


00002 I 


Wolf et al. (2007) 


2,452,522 


5455 


±0 


0001 


1949 





-0 


00013 I 


Wolf et al. (2007) 


2,452,609 


6521 


±0 


0001 


2019 








00038 I 


Baldwin & Samolyk(2004) 


2,452,684 


3143 


±0 


0003 


2079 








00012 I 


Wolf et al. (2007) 


2,452,690 


5364 


±0 


0002 


2084 








00034 I 


Baldwin & Samolyk(2004) 


2,452,731 


5986 


±0 


0007 


2117 





-0 


00185 I 


Baldwin & Samolyk(2004) 


2,452,899 


5915 


±0 


0002 


2252 








00012 I 


Wolf et al. (2007) 


2,452,964 


2991 


±0 


0001 


2304 








00001 I 


Wolf et al. (2007) 


2,452,974 


2543 


±0 


0002 


2312 








00017 I 


This paper (SOAO) 


2,452,986 


6982 


±0 


0002 


2322 








00028 I 


Baldwin & Samolyk(2004) 


2,453,019 


0518 


±0 


0001 


2348 








00000 I 


This paper (SOAO) 


2,453,055 


1393 


±0 


0002 


2377 








00046 I 


This paper (SOAO) 


2,453,062 


6057 


±0 


0002 


2383 








00058 I 


Baldwin & Samolyk(2004) 


2,453,323 


9256 


±0 


0001 


2593 








00022 I 


Baldwin & Samolyk (2007) 
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Table 3 — Continued 



HJD 




Error 


Epoch 


O-Cfuu Min 


References 


2,453,350 


0577 


±0 


0001 


2614 





0.00027 I 


This paper (SOAO) 


2,453,353 


1692 


±0 


0006 


2616 


5 


0.00081 II 


This paper (SOAO) 


2,453,360 


0128 


±0 


0001 


2622 





0.00030 I 


This paper (SOAO) 


2,453,361 


2570 


±0 


0002 


2623 





0.00011 I 


This paper (SOAO) 


2,453,387 


3896 


±0 


0005 


2644 





0.00065 I 


Hiibscher et al. (2005) 


2,453,388 


6336 


±0 


0001 


2645 





0.00027 I 


Smith & Caton (2007) 


2,453,407 


2977 


±0 


0006 


2660 





-0.00139 I 


Hubscher et al. (2005) 


2,453,407 


2987 


±0 


0035 


2660 





-0.00039 I 


Hiibscher et al. (2005) 


2,453,410 


4121 


±0 


0003 


2662 


5 


0.00205 II 


Hubscher et al. (2005) 


2,453,414 


1438 


±0 


0009 


2665 


5 


0.00060 II 


This paper (SOAO) 


2,453,425 


3416 


±0 


0001 


2674 


5 


-0.00106 II 


Wolf et al. (2007) 


2,453,641 


2436 


±0 


0003 


2848 





0.00016 I 


This paper (CbNUO) 


2,453,675 


4626 


±0 


0003 


2875 


5 


-0.00144 II 


Biro et al. (2007) 


2,453,694 


7519 


±0 


0002 


2891 





-0.00012 I 


Baldwin & Samolyk (2007) 


2,453,704 


7068 


±0 


0001 


2899 





-0.00031 I 


Nelson (2006) 


2,453,713 


4178 


±0 


0001 


2906 





-0.00001 I 


Wolf et al. (2007) 


2,453,746 


3945 


±0 


0002 


2932 


5 


0.00047 II 


Wolf et al. (2007) 


2,453,764 


4382 


±0 


0005 


2947 





0.00057 I 


Hiibscher et al. (2006) 


2,453,769 


4149 


±0 


0001 


2951 





-0.00027 I 


Wolf et al. (2007) 


2,453,793 


0585 


±0 


0001 


2970 





-0.00001 I 


This paper (CbNUO) 


2,454,054 


3794 


±0 


0001 


3180 





-0.00029 I 


Dogru et al. (2007) 


2,454,070 


5565 


±0 


0002 


3193 





-0.00022 I 


Wolf et al. (2007) 


2,454,084 


2455 


±0 


0002 


3204 





0.00052 I 


This paper (CbNUO) 


2,454,085 


4892 


±0 


0006 


3205 





-0.00017 I 


Hubscher & Walter (2007) 


2,454,115 


3547 


±0 


0010 


3229 





0.00004 I 


Hiibscher (2007) 


2,454,120 


9528 


±0 


0005 


3233 


5 


-0.00160 II 


This paper (CbNUO) 


2,454,141 


4868 


±0 


0002 


3250 





0.00002 I 


Wolf et al. (2007) 


2,454,152 


6860 


±0 


0001 


3259 





-0.00027 I 


Baldwin & Samolyk (2007) 


2,454,171 


3512 


±0 


0009 


3274 





-0.00088 I 


Hiibscher (2007) 


2,454,171 


3516 


±0 


0001 


3274 





-0.00048 I 


Wolf et al. (2007) 


2,454,176 


3298 


±0 


0001 


3278 





0.00018 I 


Wolf et al. (2007) 


2,454,186 


2843 


±0 


0002 


3286 





-0.00042 I 


Wolf et al. (2007) 


2,454,487 


4260 


±0 


0001 


3528 





-0.00045 I 


Borkovits et al. (2008) 
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Table 3 — Continued 



HJD 


Error 


Epoch 


O-Cfuii 


Min 


References 


2,454,501.11462 


±0.00007 


3539.0 


-0.00009 


I 


This paper (CbNUO) 


2,454,506.09212 


±0.00007 


3543.0 


-0.00014 


I 


This paper (CbNUO) 


2,454,510.4465 


±0.0003 


3546.5 


-0.00112 


II 


Borkovits et al. (2008) 


2,454,558.3573 


±0.0003 


3585.0 


0.00077 


I 


Brat et al. (2008) 


2,454,815.94500 


±0.00007 


3792.0 


0.00029 


I 


This paper (CbNUO) 


2,454,834.6136 


±0.0001 


3807.0 


0.00308 


I 


Brat et al. (2009) 


2,454,843.3212 


±0.0002 


3814.0 


-0.00003 


I 


Brat et al. (2009) 



Table 4. Parameters for the quadratic plus LTT ephemeris of CL Aur. 



Parameter 


Value 


Unit 


To 


2,450,097,27082±0.00046 


HJD 


P 


1.24437498±0.00000017 


d 


A 


(2.44±0.34)xl0-^° 


d 


ai2 sin ^3 


2.38±0.12 


AU 


u 


218.9±2.7 


deg 


e 


0.337±0.053 




n 


0.04556±0.00029 


deg d~^ 


T 


2,444,072±56 


HJD 


P3 


21.63±0.14 




K 


0.01330±0.00069 


d 


/(M3) 


0.0290±0.0015 


Mo 


M3 (is =90 deg)'f 


0.83 


Mq 


Ms (is=60 deg)t 


0.98 


Mo 


Ms (is=30 deg)t 


1.92 


Mo 


dP/dt 


(1.43±0.20)xl0-^ 


d yr-^ 


dM/dt 


1.30x10-^ 


Mo yr-i 



''Masses of the third body for different values of i^. 
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Table 5. Applegate-model parameters for the cool secondary of CL Aur. 



Parameter 


Value 


Unit 


AP 


1.1372 


s 


AP/P 


1.06 X 10"^ 




AQ 


8.49 X 10^° 


g cm^ 


AJ 


1.46 X 10^^ 


2 —1 

g cm s 


Is 


5.47 X 10^^ 


g cm^ 


An 


2.68 X lO"'^ 


s-i 


An/n 


4.59 X 10-3 




AE 


7.85 X lO^'i 


erg 




3.61 X 10^3 


erg 




0.926 


Le 




0.103 




ATTlrms 


±0.018 


mag 


B 


8,900 


G 



Table 6. Photometric solutions of CL Aur. 



Parameter 


Without spot 


With spot 


Primary 


Secondary 


Primary 


Secondary 


To (HJD) 


2,452,974.25441±0.00006 


2452974.25438±0.00005 


P(d) 


1.24438169±0.00000025 


1.24438175±0.00000025 


q 


0.5944(55) 


0.6023(49) 


i (deg) 


88. 


15(11) 


88.21(12) 


r(K) 


9,420 


6,315(46) 


9,420 


6,323(44) 


n 


3.569(14) 


3.053 


3.563(13) 


3.068 


A 


1.0 


0.5 


1.0 


0.5 


9 


1.0 


0.32 


1.0 


0.32 


X 


0.612 


0.478 


0.612 


0.478 


XB 


0.601(28) 


0.375(103) 


0.585(27) 


0.368(101) 


xv 


0.515(30) 


0.539(70) 


0.511(27) 


0.538(68) 


XR 


0.408(31) 


0.528(54) 


0.409(28) 


0.529(53) 


XI 


0.329(30) 


0.397(45) 


0.332(28) 


0.398(44) 


L/{Li + L2)b 


0.8629(24) 


0.1371 


0.8630(22) 


0.1370 


L/{Li + L2)v 


0.8163(24) 


0.1837 


0.8163(22) 


0.1837 


L/{Li + L2)r 


0.7683(25) 


0.2317 


0.7683(23) 


0.2317 


L/{Li + L2)i 


0.7185(26) 


0.2815 


0.7185(23) 


0.2815 


r (pole) 


0.3327(17) 


0.3135(7) 


0.3343(15) 


0.3145(7) 


r (point) 


0.3643(26) 


0.4467(29) 


0.3671(24) 


0.4481(25) 


r (side) 


0.3432(19) 


0.3276(8) 


0.3451(17) 


0.3287(7) 


r (back) 


0.3547(20) 


0.3598(8) 


0.3569(20) 


0.3609(7) 


r (volume) 


0.3439 


0.3350 


0.3458 


0.3361 


Colatitudc (dcg) 






71.87(85) 




Longitude (deg) 






2.36(29) 




Radius (deg) 






14.14(18) 




^spot /^local 






1.125(9) 




TW{0 - Cf 


0.012 


0.011 
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Table 7. Minimum timings of CL Aur determined by two different methods. 



T/" TIT" 

KvW 


W-D 


Error 


Ditterence 


t liter 


TV /T" 

Mm 


2,452,974.2543 


2,452,974.25469 


±0.00007 


-0.00039 


BVRI 


I 


2,453,019.0518 


2,453,019.05186 


±0.00010 


-0.00006 


BVRI 


I 


2,453,055.1393 


2,453,055.13934 


±0.00010 


-0.00004 


BVRI 


I 


2,453,350.0577 


2,453,350.05766 


±0.00005 


±0.00004 


BVRI 


I 


2,453,353.1692 


2,453,353.16967 


±0.00031 


-0.00047 


BVRI 


II 


2,453,360.0128 


2,453,360.01267 


±0.00009 


+0.00013 


BVRI 


I 


2,453,361.2570 


2,453,361.25700 


±0.00009 


±0.00000 


BVRI 


I 


2,453,414.1438 


2,453,414.14313 


±0.00042 


+0.00067 


BVRI 


II 



^Errors yielded by the W-D code. 
•^Differences between columns (1) and (2). 



Table 8. Astrophysical parameters for CL Aur. 



Parameter 


Primary 


Secondary 


M/Mq 


2.24 


1.35 


R/Rq 


2.58 


2.51 


log g (cgs) 


3.97 


3.77 


log p/pQ 


-0.88 


-1.07 


r(K) 


9,420 


6,323 




46.96 


9.02 


B — V (mag) 


-0.01 


+0.48 


Mboi (mag) 


+0.51 


+2.30 


BC (mag) 


-0.15 


-0.01 


My (mag) 


+0.66 


+2.31 


^v,total (mag) 




+0.44 


Distance (pc) 




1,150 



^Absolute visual magnitude from both 
components. 
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Table 9. Impact position and energy of transferred matter. 



Initial velocity (km s ^) Impact position (deg) Impact energy (Lq) 



Vx 


Vy 




Longitude 


Colatitude 




3x10-^ 








10.5 


90.0 


0.69 


3x10-^ 








10.3 


90.0 


0.69 


3x10-3 








10.5 


90.0 


0.69 


3x10-1 








10.6 


90.0 


0.70 


3 








10.5 


90.0 


0.71 


30 








11.3 


90.0 


0.73 


300 








4.9 


90.0 


1.85 


450 








3.4 


90.0 


3.24 


3x10-1 





3x10-1 


10.6 


90.0 


0.70 


3 





3 


10.6 


90.3 


0.71 


30 





30 


11.4 


85.3 


0.72 


150 





150 


12.4 


65.4 


1.26 


300 





300 


20.6 


43.6 


2.96 



